EXPLICIT BOUNDS FOR SUMS OF SQUARES 



JEREMY ROUSE 



Abstract. For an even integer k, let r2k(n) be the number of representations of 
n as a sum of 2k squares. The quantity r2k(n) is appoximated by the classical 
singular series p2k(n) x n . Deligne's bound on the Fourier coefficients of Hecke 

k — 1 

eigenforms gives that r 2 /c(n) = p2k(n) + 0(d(n)n^~). We determine the optimal 
implied constant in this estimate provided that either k/2 or n is odd. The proof 
requires a delicate positivity argument involving Petersson inner products. 



1. Introduction and Statement of Results 

In Hardy's book on Ramanujan [TT], he states the following (Chapter 9, pg. 132). 

The problem of the representations of an integer n as the sum of a 
given number k of integral squares is one of the most celebrated in the 
theory of numbers. Its history may be traced back to Diophantus, but 
begins effectively with Girard's (or Fermat's) theorem that a prime 
4m + 1 is the sum of two squares. Almost every arithmetician of note 
since Fermat has contributed to the solution of the problem, and it 
has its puzzles for us still. 

If n is a non-negative integer, let 

r s(n) = x 2 , ■ ■ ■ , x s ) eZ s : x\ + x\ H V x\ = n} 

be the number of representations of n as a sum of s squares. 

The classical work that Hardy refers to includes the formulas of Jacobi giving the 
following exact formulas. Let n be a positive integer and write n = 2 a m, where m is 
odd. Then 

, . f 8ai(m) if a = , . I 16o" 3 (m) if a = 

1 2401(777,) if a > 1, [16-^ — r ^V 3 (m) if a > 1. 

The search for higher exact formulas (each involving more complicated arithmetic 
functions) for was carried out by many mathematicians. Glaisher [9] and Rankin 
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[20] were interested in these formulas where the arithmetic functions involved were 
multiplicative. 

In a different direction, Hardy [lOj and Mordell [IS] applied the circle method to give 
an approximation 

r s (n) = p s (n) + R s (n) 

where p s (n) is the "singular series" and R s (n) is an error term. Here p s {n) can be 
expressed as a divisor sum if s is even, and p s (n) x n2~ l provided s > 4. The 
contribution R s (n) is more mysterious, and Deligne's proof of the Weil conjectures 
(see [7]) implies an estimate of the form 

(1) R s (n) = 0(d(n)ni~^) 

provided s is even. The phenomena of exact formulas for r s (n) of the form r s (n) = 
p s (n) only occurs for small s. In [21], Rankin shows that R s (n) is identically zero if 
and only if s < 8. Exact formulas of a different nature were given by Milne in [T7] 
when s = 4n 2 and s = 4n(n + 1). 

The problem we study is the implied constant in equation (CQ) above. This is a 
natural question, and in [221 El 02] , the author has studied the corresponding problem 
for powers of the A function, p-core partitions (joint work with Byungchan Kim), 
and arbitrary level 1 cusp forms (joint work with Paul Jenkins), respectively. To 
prove their now famous "290-theorem" Bhargava and Hanke [I] compute this implied 
constant for about 6000 quadratic forms in four variables and use this to determine 
precisely which integers these forms represent. 

Returning to our problem, if s — 2k and k is even, we have that 

P*(n) = ^Jr 1 ^ K-iW + (- 1 + (-l) k/2+1 )°k-i(n/2) + (-l)*/=»2*oi_i(n/4)) , 

where Bk is the fcth Bernoulli number and <jfe_i(n) is the sum of the k — 1st powers of 
the positive integer divisors of n (and is hence zero if n is not an integer). Our main 
result is the following. 

Theorem 1. Suppose that k is even. If either k/2 is odd or n is odd, then we have 

/ m ( , 2H-l) fc / 2 \ x fc -i 
|r 2fc (n) - p2k(n)\ < Uk + v _ ^ J d(n) n —. 

Remark. If 2k = 4 or 2k = 8, the right hand side is zero, and we recover the exact 
formulas of Jacobi. For arbitrary even k, we have r 2 fe(l) = 4A; and P2fe(l) = 2 ^-i) J B fc • 
Thus, the inequality above becomes an equality when n — 1. This shows that the 
implied constant 

+ (2* - l)B t 
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in (JIJ) is best possible. The error term is smaller than the main term provided n^> k 2 . 
Our approach to proving Theorem [I] is as follows. If 

l + 2^<f 2 , q = e 2wiz 



U{Z) 

n=l 

is the classical Jacobi theta function, then 

oo 



6 2k {z)=Y,r2k{n)q n 



n=0 



is a modular form of weight k on r (4). If k is even, we can decompose 

(2) 6 2k (z) = a 1 E k (z) + a 2 E k (2z) + a 3 E k (4z) + J2c l 9i(z) + J2 d ^i( 2 z) + J2 ei ^ 



where 

E,Jz) = 1 - 

B k 



oo 



E k {z) = l-—Y,° k -An)q n 

n=l 

is the classical level 1 Eisenstein series, and the gi(z) are normalized newforms of level 
1, 2, or 4. We prove the following. 



Theorem 2. Assume the notation above. Then for all i, c$ > 0. 

Theorem [2] allows us to read off 

v^, , 2H-l) fc / 2 

\ \a\=y Ci = 4k + —r { — 

\ l l L-j ^ (2 k -l)B k 

from the coefficient of q on both sides of (|2]), using that a\ = j^-- 



(2 k ~l)B k - 

To prove Theorem [2] we use properties of the Petersson inner product on Mfc(To(4)) 
(see Section [2] for precise definitions). If gi{z) is a newform of level 4, then Qi{z) is 
orthogonal to every other term in the expansion ([2]). It follows that 

(9 2k , gi ) = Ci{g h gi). 

It suffices to prove that (6 2k ,gi} > 0. This Petersson inner product consists of a 
contribution from each of the three cusps of To(4): oo, 0, and 1/2. The contribution 
from oo is 

2 ^ r 2k (n)a{n) 



y: r2k{n) k a ^ n) r u k - 2 e-du. 



4irn 



v ' n=l 

Here gi(z) = Y^Li o J { n ) < f l - Our approach is to show that the main term in the above 
sum comes from n = 1. If n is fixed, r 2k {n) is a polynomial of degree 2k in n. 
We compute these polynomials explicitly, and use this to the bound the terms when 
2 < n < 2500. Next, we use a simple induction bound on r2 k (n) to show that the 



4 



JEREMY ROUSE 



terms with 2500 < n < log (A;) are small enough. Finally, we use the exponential 
decay of f^ n u k ~ 2 e~ u du when n > J^\og(k). 

The cusp at zero behaves in an essentially identical way to the cusp at oo, and the 
contribution from the cusp at 1/2 is very small, since 9{z) vanishes there. 

Remark. This result can be thought of as a refined form of the circle method. The 
Eisenstein series is the contribution of the major arcs, while the Deligne's result, 
and the bounds we give on the constants C{ can be thought of as explicit, uniform 
minor arc estimates. Further, it is plausible that the Fourier coefficients of distinct 
newforms are independent ( an assertion that could be justified under the assumption 
of the holomorphy of certain Rankin- Selberg convolutions). Combining this with the 
recent proof of the Sato-Tate conjecture (see ^) suggests that for any e > 0, there 
are infinitely many primes p so that 



|r 2 fc(p) ~ p2k{p)\ > 



2k(-lW 2 \ „ x »_i 
4fc + T-r — k e d(p)p~ . 



Remark. The proof gives more detailed information about the constants Cj in (J2J). In 
particular, if g%{z) is a newform of level 4 and k = 2 (mod 4), then 

where a ~ 0.918. If k = (mod 4), then Cj = 0. Similar, but more complicated 
results are true for the constants Ci associated with level 1 and level 2 newforms. 



An outline of the paper is as follows. In Section[2]we give precise definitions and review 
necessary background information. In Section[3]we prove a number of auxiliary results 
that will be used in the proof of Theorem |2j In Section HI we prove Theorem [2] and 
use this to deduce Theorem [TJ Finally, in Section El we address other values of k and 
n. 

Acknowledgements. The author used Magma |5J version 2.17 for computations. 



2. Background 

In this section we give definitions and review necessary background. For N > 1, let 
Mfc(r (iV)) denote the C- vector space of modular forms of weight k on r (A) := 

lie d) ^ SL^Z;) : A|ci. Let Sf-(T (N)) denote the subspace of cusp forms. 
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a b 
c d 



If / is a modular form of weight k, and a = 
determinant, define the usual slash operator by 

' az + b 



G GL 2 (Q) and has positive 



f\a = (ad-bc) k/2 (cz + d)- k f 



cz + d 



For a positive integer d, define the operator V(d) by f(z)\V(d) = f(dz). It is well- 
known (see [12], pg. 107 for a proof) that V{d) maps M k (T (N)) to M k (T (dN)) and 
Sk(T (N)) to Sk(T (dN)). For a positive integer d, define the operator U(d) by 



no 



^a{n)q n \U{d) = ^a{dn)q n . 

n=0 n=0 

If d\N, then U(d) maps Mk(T (N)) to itself and Sk(T (N)) to itself. If p is a prime 
with p\ N, define the usual Hecke opeator T(p) by T(p) = U(p) + p k ~ l V{p). 

If /) 9 ^ ^fc(ro(AQ) and at least one of / or g is a cusp form, let 

ofcc dy 

j[x-tty)g(x-tiy)y 

'H/r (Af) 



(/^) = rqT m r / Ani // f(x + iy)g(x + iy)y 

vr[SL 2 (Z) : r (A^)J JJ m /r (N) y 



denote the usual Petersson inner product. If p \ N, then the Hecke operators T(p), 
acting on Sk(T (N)), are self-adjoint with respect to the Petersson inner product. 
Moreover, if a G GL 2 (Q) and has positive determinant, then (f\a,g\a) = (f,g)- 

Let S^ ew (T (N)) denote the orthogonal complement under this inner product of the 
space spanned by all forms 

f(z)\V(d), where f(z) G Sk(T (M)), 

and we have M\N, M < N, and d is a divisor of N/M. A newform of level N is a 
form 

oo 

f(z) = Y,a(n)q n eSr(UN)) 
n=l 

that is a simultaneous eigenform of the Hecke operators T(p), normalized so that 
o(l) = 1. We have the Deligne bound 

k — 1 

\a(n)\ < d(n)n~^~ 

where d(n) is the number of divisors of n (for a detailed proof of this inequality, see 
the new book by Brian Conrad [6]). A newform f(z) of level N is also an eigenform of 

_ -1" 



the Atkin-Lehner operator W t 



N 



N 



This operator commutes with the Hecke 



operators T(p) for primes p \ N . One has more information about the coefficient a(p) 

k 

if p\N. If N = p, then a(p) = —\p^^\ where A is the eigenvalue of / under Wn. If 
p 2 \N, then a(p) = (see [2], Theorem 3). 
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The multiplicity-one theorem states that the joint eigenspaces of all T(p) (with p \ N) 
in S^ ew (T (N)) are one-dimensional. It follows from this, and the self-adjointness of 
the Hecke opeators, that if fx and f2 are two distinct newforms, then (fx, / 2 ) = 0. It 
is known (see Section 5.11 of [8]) that the Eisenstein series Ej-(z) (and E/~(z)\V(d)) 
are orthogonal to cusp forms under the Petersson inner product. 

Finally, let rj(z) denote as usual the Dedekind eta function 

v ( z ) = q 1/24 Y[(l-q n ), q = e 2Tiz . 

n=l 

We have the following well-known identities: 

9{z) 

7] 4 (2z) 

^ + 1 ^ (stt) 

(see the exercises on page 145 of [15], solutions are on page 234). 



7] 5 (2z) 



rf (z)r] 2 (Az) 
a(2n + l)q 2 

r) 8 (4:z) 



oo 
n=0 



16 



T] A (2z) 



3. Preliminary results 

In this section we prove three lemmas that will be used in the proof of the main 
results. Our first lemma proves some simple bounds on r s (n). 

Lemma 3. (1) Suppose that n is a non-negative integer. There are non-negative 
constants Ci^ n (0 < i < n) so that 



) = $^ Ci > n ( S ) f or al1 s - °- 

i=0 



r s (n 

(2) If n is fixed, T2s a ^ is a decreasing function of s, provided 2s >n + ■ 

(3) If n is a positive integer and s > 6, then 

, , ^ 3(4.11)% . 
25^ 

Proof. We prove the first statement by strong induction on n. For n = 0, we have 
r s (0) = 1 = 1 ■ Thus, co,o = 1 and the result holds. 
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Assume the result is true for all m < n. Let t be a positive integer with t < s. Then 

Lv^J 

r t (n) - r t _i(n) = 2 r t _i(n - r 2 ) 

r=l 

Lv 7 ™] n-r 2 /t — 1 

= 2 ^2 ^2 Ci ' n -'' 

r=l i=0 



Summing both sides over alH, 1 < t < s and using that J2t=i C^ 1 ) = G+i) gi ves 

Lv^J n-r 2 , v 



r=l j=0 



= 2 C i-hn-r* 
i=l \ r =l 



Since the Cj_i jn _ r 2 are non-negative, by the induction hypothesis, it follows that their 
sum is non-negative, and this proves that the result is true for n. 

To prove the second statement, it suffices to prove that each term in the expression 

s-i / j c i,n s-i 



n 2 i=0 n 2 



is a decreasing function of s. Let /(s) = ( 2s ) • tj/ 1 s )/ 2 . Then, 

+ 1 (2s + 2)(2s+l) 



/(s) (2s + 2-i)(2s + l-i) 

1 (2s + 2)(2s+l) 
~ v 7 ^ (2s + 2 - n)(2s + 1 - n) 
1 / // 

< -F= 1 + 



n \ 2s — n 



This is a decreasing function of s, and if we take s = n + ±2 1 , then 2s — n = ^ . 
and so 

^ (i + = ^ (i + - 1)) 2 = i. 

Vn \ 2s — n/ vn 
This proves that /(s + 1) < /(s), as desired. 

We prove the third statement by induction on s. Our base case is s = 6 and in this 
case, we use the exact formula 

r 6 (n) = ]T d 2 (-4 X _i(d) + 16x_i(n/d)) , 
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where 

{1 if n = 1 (mod 4) 

—1 if n = 3 (mod 4) 

if n is even. 

We rewrite this as 

If ra is even, then r 6 (n)/n 2 < 8((2) < 13.2. On the other hand if n is odd, then 

I6x-i(n/d) —Ax~i(d) is negative if n/d = 3 (mod 4) and 16x_i(n/d) — 4x_i(cf) < 20 
iin/d= 1 (mod 4). Thus, 



rein) 

n 2 d 2 -= - ^ (4n + l) 2 ' 



1 ^ 1 
< 20 > — < 20 > 

d\n 

d=l (mod 4) 



One can show that the right hand side above is about 21.4966613 < We denote 



300 

by C s a constant so that r s (n) < C s (n + and we take C 6 = This proves 

the base case. 

Assume now that s > 6. We have 

r s +i(rc) = r s (n) + 2 ^ r s (n - m 2 ) 

m=l 

< C s (n + s)^ 1 + 2C S ^(n + s-m 2 )^ 1 

m=l 



< CJra + s)^ 1 + 2C S / (n + s + 1 - a; 2 )^ 1 dx 



We have 



Thus 



< C s {n + s)^" 1 + 2C s (n + s + l)^" 1 / (1 - w 2 )^ 1 du. 

io 



(1 - u 2 )^ 1 = e (f-l)Ml- 2 ) < e -(s/2-lW_ 



2j\l- u 2 )^ 1 <*« < 2 jH e-^ 2 -^ 2 du = J^O, 
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and 



r s+ i(n) < C s {n + s)^ 1 + C s {n + s + l)^" - 1 



2tt 
s-2 



< CJn + s + 1 



8 + 1 



2tt (n + s)^- 1 
+ 



s-2 ( n + s + 1)^-1 



Note that the second term inside the brackets above is a decreasing function of n and 
is relevant only for n > 1. It follows that 



r s+ i(n) < C s {n + s + l)^" 1 



C, • 4.11 



s-2 V s + 2 



s/2-l' 



< * ' (n + s + 1)—- 1 . 
Hence, we may take C s+ \ = -^=C S and so 



_ 6449 4.11 s - 6 3(4.11) J 



300 v /iI76! ~ 25Vs! ' 



□ 



Next, we use Deligne's bound on the Fourier coefficients of a newform to bound its 
value. 

Lemma 4. Suppose that k > 7, y > and g(z) = Y^=i a(w)<Z n with \a(n)\ < 

k — 1 

d(n)n^2~ . Then 



\g(x + iy)\< ^rr^r ( " 



{2ixy) "t 1 



log ( )+7 + l 



where 7 is Euler's constant. 



k — 1 

Proof. Since the nth Fourier coefficient of g(z) is bounded by d{n)n~^~ , we have that 

00 

|s-(:c + iy)| < ^d(? 



fc-i 



n n 2 e 



-2irny 



n=l 
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If D(x) = J2 n <x d(n), then D(x) < x\og{x) + 7X + 1 < x\og{x) + (7 + l)x. By partial 
summation, we have 



n=l 



n)n 2 e 



-2irny 



D(x) 



fc-i I k — 1 \ fc-3 
2nyx 2 — ( — - — ) x 2 



e~ 27TXy dx 



POO 

< 2iry / (log(x) + (7 + l))x^ e~ 2nxy dx. 

4:7ry 

Now, we set u = 2nxy, du = 2ny dx. We get 
2ny 



k + l 

log I T^r.) +(7+1)1 I — 1 2 



2ny 



du 



2wy 



1 /"°° k 1 

= 7T+T / (log(«) - log(27ry) + 7 + 1) w^e"" du. 

{2iry)— 

Since y > log(27n/) > and so we neglect the term involving it. We get 
1 



(27n/) ^ 



00 /-0O k 1 

log(w)w~e~ u + (7 + 1) / u~e~ u du 
fe-i /fc-i 

2~ 17 



If we extend the integrals down to zero, then the negative contribution of 



C log(w)w 2~e w du is cancelled by that of [1.5, 2] for k > 7. Thus, we get the bound 



(27TJ/) ^ 



r 



fe + i 



where V(z) = J °° \og(u)u z l e u du and ip(z) = ip(z) = The formula (see 

equation 6.3.21 on page 258 of pQ) 

2t dt 



1 r 

il)(z) = log(s) - — - / 
zz Jo 



! (z 2 + t 2 )(e 2 ^ -1) 
shows that ip( z ) < l°g(+)- Thus, we obtain the bound 



I v (k + l 



(27T2/)^ 



, fc + 1 , 

log — — +7 + 1 



□ 



Finally, we will need to understand Petersson inner products of newforms / with their 
images under V(d). This is the subject of the next result. 

Lemma 5. Suppose that f(z) = J2^ =l a(n)q n e S£ ew (T (N)) is a newform. If p\ N, 
then 
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Note that the assumption that / has trivial character implies that the Fourier coeffi- 
cients of / are real. This fact will be used frequently in what follows. 

Proof. Rankin proved in [19] that if / = ^2 a(n)q n and g = b(n)q n are cusp forms 
of weight k, then 



n<x 



We will use this formula to prove the results above. We start by letting c 
and p be a prime number with p \ N . Then, 

</^iv-( P )>=iimi.Iy:«M 

x-^co c x z — ' n 

n<x 

. 1 1 a(pn)a(n) 

- hm - - — y 

x—>no r T ' J 



(fc-1)! 



x-»-oo c x ^— ' (pn) k ~ 1 

111 a(pn)a(n) 



lim ----- 



x— >oo C V K — ' Tl k 1 

— p 



= ^</,/I^(p)>- 

Now, a(p)/ = f\T(p) = f\U(p) +p k ~ 1 f\V(p). It follows that 

a(p)(f,f) = (fJ\T(p)} = (fJ\U(p))+p k - l {fJ\V(p)) 
= p k (fJ'\V(p))+p k - 1 (f,f\V(p)) 
= p k ~ 1 (p+l)(f,f\V(p)). 



Thus, 



4. Proof of Theorem CD and Theorem [2] 



□ 



In this section, we will prove the main results. We will first prove Theorem [2] and 
then deduce Theorem [1] from it. 

Proof of Theorem^ First, for each newform g of level 1, 2, or 4, we will find a form 
g with the property that the coefficient of g in the representation of 9 2k is positive if 
and only if (6 2k ,g) > 0. Each g will be an eigenform of T p for all odd primes p, and 
will also be an eigenform of W4 with the same eigenvalue as that of 9 2k . 
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Recall the decomposition 



6 2k (z) = a x E k (z) + a 2 E k {2z) + a 3 E k (4z) + ^ ag^z) + ^ d igi {2z) + ^ e^(4z), 



where the ^ are newforms of level 1, 2, or 4, and the q, di, el. If is an eigenspace 
for all T n (with n odd), then V is also stable under W4. Since 2fc |W r 4 =(-l)f0 2fc , it 
follows that the projection of # 2fc onto V must also have eigenvalue (—1)2 under W4. 

If V is an eigenspace coming from a newform gi of level 4, then dim V = 1. If q 7^ 0, 
then ft|W4 = (—1)2. In this case, we have (9 2k ,g) = (ci9u9i) — c i(9ii9i) an d thus 
q > if and only if (0 2k ,gi) > 0, and so we set g^ = g^ Part (i) of Theorem 7 of 
[2 J shows that for any newform of level 4, gi\W\ = —1, and hence q = if k = 
(mod 4). 

If is an eigenspace coming from a newform ^ of level 2, then dimV^ = 2. This 
vector space decomposes into one-dimensional plus and minus eigenspaces under the 
action of W4. It follows that the projection of 9 2k onto V is c*(c/i + (— 2) 2 A<?j| V(2)), 
where A is the Atkin-Lehner eigenvalue of g^ Thus, we set (ji = gi + (— 2) 2 Xgi\V(2). 
This form will be orthogonal to any element in the opposite W4 eigenspace, since W4 
is an isometry with respect to the Petersson inner product. It follows that q > if 
and only if (9 2k ,(ji) > 0. 

If V is an eigenspace coming from a newform g^ of level 1 , then dim V = 3 and 



We have that V = V + © V~ , where V + and V~ are the plus and minus eigenspaces 
for W4. Then dimV + = 2 and it is spanned by g { + 2*p i |V(4) and gi\V(2). Also 
dimV~ = 1 and it is spanned by g^ — 2 k gi\V(4). If k = (mod 4), then the Atkin- 
Lehner sign is +1. If k = 2 (mod 4), the Atkin-Lehner sign is —1. 

When k = 2 (mod 4), we set gi = gi — 2 k gi\V{4). This form satisfies #i|W4 = — <?j, 
and is again orthogonal to the form spanning the plus eigenspace for W4. 

When k = (mod 4), we set = |a(2)^|y(2) + 2 fc ^|\/(4). This form satisfies 
c/i\W4 = §i, and is hence orthogonal to gi — 2 k gi\V(A). By Lemma [5] it is orthogonal 
to gi \V{2). 



9i\W 4 
9i \V(2)\W4 
9i\V(4)\W4 



2^1^(4) 

9i\V(2) 

2~ k gi . 
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We have 



l1X J Ji/rn(4) 



dx dy 



H/r (4) 

1/2 poo 



— V / / {0 2k \kl 3 ) (x + iy)g i \ klj (x + iy)y k - 2 dydx. 



Here, the matrices 





1 


0" 




"0 


-1 




"0 


-f 




"0 


-f 




"0 


-1 




1 


0" 


7i = 





1 


,72 = 


1 





,73 = 


1 


1 


,74 = 


1 


2 


,75 = 


1 


3 


,76 = 


2 


1 



are a set of representatives for the right cosets of r (4) in SL 2 (Z). 

Term 1: This is the contribution from the cusp at infinity. In particular, it is the 
j = 1 term in the above sum. We split this term into two parts: {x + iy : —1/2 < 
x < 1/2, y > 1}, and {x + iy: -1/2 < x < 1/2, Vl - x 2 < y < 1}. 

Write 



9i(z) = ^a(n)q r ' 



n=l 



Applying the Deligne bound to each of the various possible forms of g~i, we see that 
in all cases \a(n)\ < ^-d(n)n~ . 

The first part is 



— J J iJ2 r 2k {m)e-^ my e 2 ~\ f ^ a(n)e- 2my e~ 2mnx J y k ~ 2 dx dy. 



Since the Fourier series representations converge uniformly on compact subsets of 
these regions, we can invert the summations and the integrals and obtain 



1 00 00 p i -/' z r°° 

™— n — 1 J — 1/2J1 



m=0 n=l 



•1/2 /-oo 
-1/2 



The integral over — l/2<x<l/2is zero unless m = n, in which case it is 1. We set 
u = 4nny, du = Aim dy and this gives 
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We now split this sum into several ranges. The main contribution comes from n = 1. 
We have a(l) = 1 and r2fc(l) = Ak. This term is 



8A; 



u k - 2 e~ u du 



In 



> 



8k 

8k 
(4^> 



u 



k - 2 e~ u du - 



An 



u k ~ 2 e- u du 



[(Jfe-2)! - (4vr) fc - 1 e- 4 " 



for k > 15, since if k > Air + 2, u 2 e u is increasing on [0, Air]. 

The second range is 2 < n < 2500. Here we explicitly compute the polynomials r2k{n) 
(using the algorithm in the proof of part 1 of Lemma [3]) . Part 2 of Lemma [3] shows 
that r2 jff is a decreasing function of k, provided k > 1456. 

The third range is 2500 < n < ^- log(2/c). In this range, we use the bound from part 
3 of Lemma [3j the Deligne bound, d(n) < 1\fn 1 and we obtain that 



r 2k (n a(n 



n 



k-l 



< 



< 



34(4.11) 



2 k 



25^/(2k)\ 

34 (4.11) 2fc 
25 TPfel V 2tt 



(n + 2k) 

n 2 



fe-i 



log(2fc) ■ v 7 ^ + 



2A; 



k-l 



The function /(x) = (x + ^) 
ve tha 

bound 



2k\ k ~ 1 



is decreasing for x < \2k and increasing after that. 
We have that /(50) = /(§) and § > log(2fc) if A; > 724. Thus, we have the 



68 (4.11 



fe 3/2 



25 (47r) fc v ^2fc)! (2tt) 3 / 2 



log 3/2 (2A;) • ( 50 + 



2A; 
50 



fe— l 



(k -2)!, 



valid provided k > 724. For k < 724, we use the larger of the values of / at x = 50 
and x = J^\og(2k). 

The fourth and final range is n > ^log(2A;). In this range we use the decay of the 
integral J e 



Ann 



U 



" 2 e u du. We have that u > 2k log(2fc) and so u k 2 e u < e u l 2 and so 

the integral is bounded by 2e~ 2wn . Bounding a(n) and r2/c(n) as before, we have that 
the contribution from this range is at most 



34 



3(4vr) 



2nl 



3 (4.11 



,2k 



n=Alog(2fc) 



n*- 1 125 v /p)T 



(n + 2A;) fe - 1 ■ 2e~ 2 ™ 
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We write 



(n+2k) 



r-r— as (1 + f ) k \ If k > 40, 1 + f < 3.87 and we get 



136 (4.11) 2fc (3.87)^ 1 



25{Att) 



oo 

E 

n=£log(2fc) 



n 2 e 



-2irn 



If a n — n^e 2?m , then we have 



*n+l 



< 1 + 



n 



-2tt 



< e ^' 2n < e < e 



27T+, 



-5.6 



Thus, we get the bound 



136 (4.11H3.87)'- f}_ log{2k)] - {2k) -* 



25(4vr) fc 
valid if k > 40. 

The second part of the contribution from the cusp at infinity is 

i r 1/2 r 1 



l-e- 



2tt 



9 2k (x + iy)g(x + iy)y k 2 dydx. 

-1/2 J^l-x 2 

In this region we use Lemma H] to bound g(x + iy), and we use that 



\6{z)\ < l + 2]TV 2 ™ 2y < 1.008667 



n=l 



for y > y/3/2. This gives the bound 

r (^) [log (^) + 7 + 1] (1.008667) 2fc f 1 ' 2 f 1 



(27TJ 2 



y~ dydx. 



J -1/2 Jy/T^x 7 

The double integral above is less than or equal to f^ 2 Jq y^ dydx = Hence, 



we get the bound 



valid for k > 7. 



34r (*±1) [log (*±1) +7 + 1] (1.008667) 2fc 



3(fc-3)(27r 



fc+3 



Term 2: This is the contribution of the cusp at zero, and in particular the contribu- 
tions from the terms involving 72, 73, 74, and 75. We have 

9 2k \W 4 = (-l)^ 2fc and gi\W 4 = 
Translating this into Fourier expansions gives 



i2k\ 



-1 

1 



2 fc V4 





"0 -1" 


) ) &l 


1 



4 



2 fc & 
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Thus, the contribution from these four terms is 
3 r i/2 r °° 

E ' 



1 



(2?r) • 4 fe 



3=0 



e 



2k 



x + j + iy 



9i 



-1/2 JVl^x 1 

We set u = x/4 and v = y/4 in the integrand and obtain 

3 



' y k 2 dy dx. 



-Erf ^L +iv+ i) 



Pi ( w + if + ^ ) v k 2 dv du. 



We break this into two terms. The first term consists of those pieces with v < 1. 
The smallest value v takes on this piece is \/3/8 and since a/3/8 > ^ we may use 
Lemma H] to bound the contribution. This yields 



\gi(u + iv)\ < 



17 



. r( *±l 

3- (2th;) 4 * 1 V 2 



log 



k + 1 



+ 7 + 1 



We also have 



+ < 1 + 2^ e" 2 ™ 21 ' < 1.52182 

n=l 

for v > v3/8. The contribution of these terms is therefore bounded by 
17- (1.52182) 2fc 



3- (2tt)^ 



-r 



k + l 



log 



*; + i 



+ 7 + 1 



3 

E 

3=0 



1/8 



\/l-16n 2 

4 



f 2 (if 



The sum of double integrals is bounded by f_y 2 f * 2 ° dv = and we get the 



bound 



2k 



34- (1.52182) 
3- (2vr)^ . (jfc-3)' 



k + l 



l0£ 



fc + 1 



+ 7 + 1 



on the part where v < 1, valid for k > 7. 

The second term consists of those pieces with u > 1. This gives 



1 

2^ 



1/2 /-co 



1/2 Jl 



8 (u + iv)gi(u + iv)v dv du. 



This is exactly the same as the contribution of the first part of the cusp at infinity! 

Term 3: This is the contribution of the cusp at 1/2 corresponding to the matrix 7@. 
We must understand the Fourier expansion of gi\j6. Since 7 6 G r (2), terms of level 
1 or level 2 are not affected. 

If g is a newform of level 4, then since 76 is not in To(4), we have that g h-> g + g|76 is 
the trace map from 5^ (To (4)) to Sk{To(2)). Since newforms are in the kernel of the 
trace map (by Theorem 4 of [IE]), it follows that g + gl^e = and so g|76 = —g. 
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g\V(2)\ l6 = 2- k / 2 



9\ 



2 
1 



1 

2 1 



2 -ft/2 







-1 
2 



1 
-2 1 



1 1/2 
1 



The first matrix is the Atkin-Lehner involution of level 2, of which g is an eigenform. 
The second matrix is in To (2) and the third matrix does not affect the size of the 
Fourier coefficients at infinity. It follows that the nth Fourier coefficient of <?|V(2)|7 6 
is bounded by 2~ k l 2 d{n)rr~^ . 

If g is a newform of level 1, we have 



g\V{A)\ 1(i = 2~ k g 



"4 0" 




"l 0" 




"2 l" 




"2 -1" 


1 




2 1 


= 2~ k g\ 


1 1 




2 



2-^ -1/2). 



Thus, the nth Fourier coefficient of #1^(4) I73 is bounded by 2 d(n)n 2 . It follows 
that for any c/i, the nth coefficient of chI'Jg is bounded by ^-d( 

Now, 6 2k \ 1& = 2 2k ^0^. The form F(z 



k — 1 

n)n~ 



r]{Az) s 



v (2zY G ^2(r (4)) and satisfies 
F(z) = 



n odd 



Thus, for y > y/3/2, \F(z)\ < e~ 2 ^ (E„ odd cr(n)e- 2 ^ n -^ y ) < 1.0001e~ 2 ^ and so 



ft/2 p-krty 



\0 2k \l$\ < 2 2fe (1.0001) fe/2 e 
for y > y/3/2. The contribution of the cusp at 1/2 is therefore 

-1/2 POO 



1 

2^ 

By Lemma HI we have 



1/2 JVT^? 



6 |7 6 (x + iy)gi\^(x + iy)y dydx. 



\9i\l&{x + iy)\ < 



This gives the bound 



14 



(2vr)^ 



r 



log 



jfc + 1 



+ 7 + 1 



14-2 2fc - (l.OOOl)*/ 2 fk+1 
3(2vr)^ 



log 



k + 1 



+ 7 + 1 



y 2 



The integral above is — tt^-T (J*yj, and so the bound on this term is 



U-2 2k ■ (1.0001) fc / 2 r (k I 



3(27r)^(/br)^ 
and is valid for k > 7. 



r 



jfe-3 



log 



k+ 1 



+ 7 + 1 
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After dividing each term above by ^ 4 Jy!' , the main term is increasing linearly (it 
is about 16k), and each other term decreases exponentially. The most troublesome 
term is the term from the third range of values of n from the cusp at infinity, and 
(after dividing by felr) is asymptotic to Cic^k 1 ^ 4 ln(fc) 3//2 , where c 2 ~ 0.918, but 
C\ ~ 1.69 • 10 543 . This term is larger than the main term only when k > 14000. 

For this reason, we must explicitly calculate our bounds for k < 14000. In this 
range, we refine our estimate of the troublesome term by using the exact values of the 
incomplete T-function u k ~ 2 e~ u du. Also, for k < 2550, we compute the values of 
r 2 fc(n) explicitly for 2 < n < 2500 and use these in our bounds. For k > 2552, we use 
part 2 of Lemma |3j 

Finally, for k < 194, our numerical bounds are not sufficient and we use Magma 
to explicitly compute the decomposition of 9 2k as in equation (j2]) and find that the 
constants q are non-negative. □ 

Proof of Theorem^ First, assume that n is odd. Considering the coefficient of q on 
both sides of ()2]), we obtain 

2k(-l) k ' 2 ^ 

By Theorem [2], we have 

2A;(-l) fc / 2 
(2 k -l)B k ' 



Deligne's bound on the nth coefficient of g%{z) is bounded by d(n)n^~ . Plugging this 
bound into the decomposition and using the fact that the coefficients of q n in gi{2z) 
and gi(iz) are zero if n is odd gives the desired bound on the cusp form contribution 
to d 2k (z). 

Now, suppose that k/2 is odd and n is even. Then k = 2 (mod 4). We represent the 
decomposition of the cusp form part of 6 2k (z) as 

C{z) = 5> (/<(*) - 2 k f i (z)\V(A)) (W) - 2h l9t (z)\V(2)) 

i i % 

Here, the fi(z), g%{z) and hi(z) are the newforms of levels 1, 2, and 4 respectively, 
and Aj is the Atkin-Lehner eigenvalue of gi(z). One can see that the nth coefficients 
of fi(z) - 2 k fi(z)\V(A) and gi (z) - 2? X igi (z)\V(2) are bounded by 3d(n)n 1 ^ 1 . Thus, 
for even n, we obtain the bound 



3r\ + 3,Si J d{n)n 



2 



We will show that £\ 3rj + 3s; < 4/c 
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2k 



(2*-l)B fc ■ 



To compute the constant £V 3rj + 3sj, we will compute the trace of C(z) to 5 , fc(r (2)), 

' . Straight-forward, but somewhat lengthy 



2 1 



given by Tr(C7) := C{z) + C{z) 
computations show that 

Tr(/,(z) - 2 k f l {z)\V{A)) = - 2o i (2)/ i (*)|V(2) 

Tr(g i (z)-2%\ i g i (z)\V(2) = 3g i (z) 
Tr(hi(z)) = 0. 

It follows from these formulas that + 3s« is the coefficient of q in Tr(C). We 

have that 



Tr(C) = Tr(£ 2fc ) - ^J—Tr(E k (z)) - ^—^{E^z)) 

= {f + 4k p^) + 2*3I £ *W - (C 1 + 2 1 -*)^(z)|V(2) - 2-*£ fc (* 

Taking the coefficient of q on both sides of the preceding equation gives 

> 3r, + 3sj = 4k — -t— : r-=— < Ak — -7—, r-=— 

Y (2 fc -l)5 fc (2 fc -l)B* 

since k = 2 (mod 4) and hence > 0. This proves Theorem [1] in the case that 
k = 2 (mod 4) and n is even. □ 

5. Final remarks 

It is natural to consider if Theorem [1] is true in other cases. When k = (mod 4) 
and n is even, the main issue is that if gi is a level 1 eigenform and 

. 00 
h = 9i~ -a{2) 9i \V{2) + 2 k gi \V{A) = Vc(n) 9 " 

6 n=l 

k — 1 

then the best possible bound on the Fourier coefficients of <?j is \c(n)\ < 2d(n)n~. 
In order for this bound to come close to being achieved, it is necessary for |a(2)|, the 

k — 1 

absolute value of the second coefficient in g^ to be close to 2 • 2~^~ . Serre proved 
in 1997 (see [22]) that if p is a fixed prime, the pth coefficients of newforms become 
equidistributed (along any sequence of weights and levels whose sum tends to infinity, 
where the levels are not multiples of p). It follows from this that there will be level 

k — 1 

1 eigenforms with |a(2)| arbitrarily close to 2 • 2^~, but also that there will be few 
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such forms. One approach to extending Theorem [T] to the case when k = (mod 4) 
is to use the equidistribution of the numbers |a(2)|. 

It is also natural to consider the problem of deriving a sharp bound in the case that 
k is odd. In the case when k is even, the contribution from the cusp at zero is (up to 
a fairly small error) the same as the contribution at the cusp at infinity, since both 
8 2k and the newforms are eigenforms of the Atkin-Lehner involution W4. However, 
when k is odd, the newforms are not eigenforms of W4 any longer. This means that 
the contribution of the cusp at zero is (up to some small error) the contribution of 
the cusp at infinity times some complex number A of absolute value 1. This complex 
number is related to the coefficient of q A of the relevant eigenform A similar result 
could be proven provided one could rule out the possibility that A is close to — 1. 
In fact, the analogue of Theorem [2] is false for k = 17, although this seems to be a 
consequence of the smallness of the weight, rather than a value of A too close to —1. 

For half-integral values of k (corresponding to representations of n as the sum of an 
odd number of squares), the question is still interesting. In this case, the coefficients 
of the cusp forms involve square roots of central critical L-values of quadratic twists 
of forms of level 1 and level 2. The analogue of Deligne's theorem in this case would 
be optimal subconvexity bounds on these L-values, currently attainable only under 
the assumption of the generalized Riemann hypothesis. 
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